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Abstract 

We investigate the possibility of testing of the Einstein Equivalence Princi- 
ple (EEP) using measurements of anomalous magnetic moments of elementary 
particles. We compute the one loop correction for the g — 2 anomaly within the 
class of non metric theories of gravity described by the THefi formalism. We 
find several novel mechanisms for breaking the EEP whose origin is due purely 
to radiative corrections. We discuss the possibilities of setting new empirical 
constraints on these effects. 
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I Introduction and Summary 



Metric theories of gravity offer the singular beauty of endowing spacetime with a sym- 
metric, second-rank tensor field g^ v that couples universally to all non-gravitational 
fields. This unique operational geometry is dependent upon the validity of the Ein- 
stein Equivalence Principle (EEP), which states that the outcomes of nongravitational 
test experiments performed within a local, freely falling frame are independent of the 
frame's location (local position invariance, LPI) and velocity (local Lorentz invariance, 
LLI) in a background gravitational field. Non metric theories break this universality 
by coupling additional gravitational fields to matter, and so violate either LPI, LLI or 
both. Limits on LPI or LLI are imposed by gravitational redshift and atomic physics 
experiments respectively, by comparing atomic energy transitions that are sensitive 
to these symmetries. Laser experiments have set stringent limits on violations of LLI 
(to a precision of ~ 10~ 22 ) [TJ, while the next generation of gravitational redshift 
experiments could reach a precision up to 1CT 9 M. 

These experiments probe transitions that are predominantly sensitive to nuclear 
electrostatic energy, although violations of WEP/EEP due to other forms of energy 
(virtually all of which are associated with baryonic matter) have also been estimated 
0. However there exist many other physical systems, dominated by primarily non- 
baryonic energies, for which the validity of the EEP is comparatively less well under- 
stood Q]. Included in the list of such systems are photons of differing polarization 
0, antimatter systems 0, neutrinos mesons ||, massive leptons, hypothesized 
dark matter, second and third generation matter, and quantum vacuum energies. 

In order to establish the universal behavior of gravity, it is important to empiri- 
cally confront the EEP with as diverse a range of non-gravitational interactions as is 
possible. Amongst the most intriguing of the list of physical systems noted above are 
those for which potential violations of the EEP can arise in quantities dependent upon 
vacuum energy shifts, which are peculiarly quantum-mechanical in origin (i.e. do not 
have a classical or semi-classical description). The empirical validity of the EEP 
in physical systems where radiative corrections are non-negligible remains an open 
question. Such systems provide an interesting empirical regime for a confrontation 
between gravitation and quantum mechanics. 

Perhaps the most useful framework for further such testing of the EEP is provided 
by quantum electrodynamics (QED). We considered this approach in a previous paper 
@ by analyzing the behavior of Lamb shift transition energies within the context of a 
wide class of nonmetric theories of gravity. The Lamb shift, along with the anomalous 
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magnetic moments (g—2 factors) of elementary particles, presents the most compelling 
evidence in support of QED. It is therefore natural to extend the nonmetric analysis 
to the g — 2 anomaly. 

We consider in this paper the possibility of using measurements of anomalous mag- 
netic moments of elementary particles as a possible test of the EEP. The high precision 
attained in g — 2 experiments motivated the early work of Newman et al. (K| to use 
such experiments to set new bounds on the validity of special relativity. Similarly, 
we expect (and shall subsequently demonstrate) that such experiments could impose 
stringent and qualitatively new limits on the parameter space of nonmetric theories. 
Our interpretation of those experiments is substantially different from theirs as we 
explicitly include violations of the EEP in the computation of g — 2. These effects 
were not considered in ref. ||10|| , which assumed that violations of special relativity 
arose only from the dynamical equation governing the motion of the fermion. 

We consider the class of non-metric theories described by the Tife/i formalism 
||11|| , following the approach given in ref. |J in developing gravitationally modified 
(GM) QED. This formalism encompasses a wide class of nonmetric theories of gravity, 
and deals with the dynamics of charged particles and electromagnetic fields in a 
static, spherically symmetric gravitational field. It assumes that the (classical) non- 
gravitational laws of physics can be derived from an action: 

Sng = -£>« f dt(T-Hv 2 a ) 1/2 + Y.^ f dtv»A^ a ) 

+±Jd 4 x(eE 2 -B 2 /^, (1) 

where m a , e a , and x%(t) are the rest mass, charge, and world line of particle a, x° = t, 
v% = dx£/dt, E = — VA — dA/dt, B = V x A. The metric is assumed to be 

ds 2 = T(r)dt 2 - H(r) (dr 2 + r 2 dtt 2 ) (2) 

where T, H, e, and /x are arbitrary functions of the (background) Newtonian gravi- 
tational potential U = GM/r, which approaches unity as U — ► 0. For an arbitrary 
non-metric theory, these functions will depend upon the type of matter, i.e. the 
species of particle or field coupling to gravity. The functions e and /i parameterize 
the 'photon metric', whereas T and H parameterize the 'particle metric' in the static, 
spherically symmetric case. Although we shall generically employ the notation T and 
H throughout this paper, it should be kept in mind that these functions shall in gen- 
eral have one set of values for electrons, another set for muons, another for protons, 
etc.. Universality of gravitational coupling in the particle sector implies that the T 
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and H functions are species independent. It is an empirical question as to whether or 
not such universality holds for all particle species. The stringent limits on universality 
violation set by previous experiments |l] have only been with regards to the relative 
gravitational couplings in the baryon/photon sector of the standard model. For the 
leptonic sector relevant to our considerations, relatively little is known jp. 

In order to study atomic systems, this classical action can be generalized for 
quantum mechanical systems ||12||. 



S = J d 4 xtfj(i ^ + e4-m)ip + -J d 4 x{E 2 -c 2 B 2 ), (3) 

where local natural units are used, A = r y fl A^, and c 2 = H /T eo[io with the sub- 
script "0" denoting the functions evaluated at X = 0, the origin of the local frame 
of reference. This action emerges upon replacing the point-particle part of the action 
in (|l]) with the Dirac Lagrangian, expanding the THefi parameters about the origin, 
neglecting their spatial variation over atomic distance scales, and rescaling coordi- 
nates and fields. These operations have the effect of treating spatial and temporal 
derivatives and gauge couplings in the fermionic sector on the same footing, and so 
we need not add additional terms of the form ip(i^f°(dQ + eAo))i/) times an arbitrary 
constant, as such terms can be re-absorbed into a definition of the parameters which 
appear in (||) above, as noted below. 

The action (§) refers to the preferred frame, as defined by the rest frame of the 
external gravitational field U . In order to analyze effects in systems moving with 
velocity u with respect to that frame we need to transform the fields and coordinates 
in (IID via the corresponding Lorentz transformations. This gives 



S = j ' d A xi){i $ + e - m)ilj + J ' d 4 x 



1 



+ - / d*x (E - B 2 ) (4) 



2 

+ £7 2 (u 2 E 2 - (u ■ E) 2 + B 2 - (u ■ B) 2 + 2u ■ (E x B 

where is the electromagnetic 4-current associated with some external source, 7 2 = 
(1 — -u 2 ) -1 , and £ = 1 — c 2 is a dimensionless (and species- dependent) parameter that 
measures the degree to which LPI/LLI is broken. This parameter scales with the 
magnitude of the dimensionless Newtonian potential, which is expected to be much 
smaller than unity for actual experiments. We are therefore able to compute effects 
of the terms in eq. ([|) that break local Lorentz invariance via a perturbative analysis 
about the familiar and well-behaved c — > 1 or £ — > limit. 
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We therefore consider gravitationally modified Quantum Electrodynamics (GMQED) 
based on the action (|j). The fermion sector and the interaction term do not change 
with respect to the metric case, and so neither do the fermion propagator and the 
vertex rule. All the non-metric effects are accounted for in the pure electromagnetic 
sector of the action, which in turns modifies the photon propagator. After a proper 
choice of the gauge fixing term involved in the quantization of the photon field, it 
can be shown || that the photon propagator takes the form up to first order in £ (in 
momentum space): 



^ = -(1 + 0^ + 4 



(5) 



where r]^ is the Minkowski tensor with a signature (H ); 7 2 = 1/(1 — u 2 ) and 

/3 M = (1, u); henceforth j3 2 = 1 — u 2 . 

Therefore eq. (|5|), along with the unmodified fermion propagator Sf{p) and vertex 
rule, may be used as the basis of the Feynman rules of GMQED. The radiative 
corrections affecting these quantities are defined in terms of the photon self energy 
IP" (A;), fermion self energy S(p), and vertex function T M respectively. These insertions 
involve the calculation of loop integrals as given by the Feynman rules up to a given 
order. 

The addition of more parameters to the theory also entails new renormaliza- 
tions beyond those of the wavefunctions, charge and mass of the fermion. The 
Tife/i parameters appear as functions of Cg = T /H and c 2 = l//ioeo, and must 
then be correspondly redefined. In units where Co = 1 (c* = c), EEP- violating correc- 
tions only appear in the electromagnetic sector of the action (as terms proportional 
to £). However a more general choice is Cq ^ 1, for which the particle sector of the 
Lagrangian density is of the form 

C D = y-m)il> + £ ?(po - V&hV (6) 

with £ = 1 — c o ■ m the moving frame this is (up to a constant) 

C' D = ^-f-m)^ (7) 

+ e 7 2 W -p-P- v) h . 

From (0) we see that quantum corrections of the form 

5C D = ^ 1] (3 ■ p - 6^ 2) (3 ■ V) 
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can still be expected. It is straightforward to show that gauge invariance will guar- 
antee 5$> = 5^ 2) = 5£ Q . 

Hence, in order to renormalize the mass and the THefi parameters, we have to 
include counterterms of the form 

5m + 5£ (S(P-p-P-V), (9) 

which consequently participate in the redefinition of the fermion self energy and vertex 
function. Additionally, the electromagnetic sector will induce quantum fluctuations 
of the form: 

5C EM = 5^{(k 2 - 7 2 (/3 • k) 2 ) Vllu - 7 2 /yU 2 }^ (10) 

corresponding to the renormalization of the THe/u, parameters, or equivalently £ 
1 — H /T fi eo, entailing a renormalization of the photon self energy. 

This summarizes the procedure for performing perturbative calculations in GMQED 
as employed previously in ref. |§ in computing the Lamb shift. This calculation was 
complicated by the boundedness of the electron to the hydrogenic atom under con- 
sideration. Here, in the g — 2 case, we have to deal only with free leptonic states. In 
Sec. II we evaluate the (one loop) radiative corrections to the elastic scattering of a 
free electron by an external electromagnetic field. Sec. Ill relates the scattering am- 
plitude with experimental observables describable in terms of the g-2 anomaly. This 
requires a derivation of the relativistic equation of motion for the electron spin in the 
presence of a magnetic field. This follows from a classical treatment for the electron 
and magnetic field. Quantum effects are introduced in the modified Hamiltonian only 
(as are non metric effects). The connection with possible LLI/LPI violating experi- 
ments is presented at the end of this section. Discussion of the results and general 
comments are given in Sec. IV. More details about the loop calculation are shown in 
the Appendix. 



II (GM) Free Scattering 

We shall consider the lowest order radiative correction to the elastic scattering of elec- 
trons by a static external field N 1 . These one loop contributions can be summarized 
in terms of the Feynman diagrams illustrated in Fig. 1. 

The Feynman amplitudes for the diagrams follow from the Feynman rules giving 
the result [|L3[| : 

A"(p',p) = u{p) {r^ + + L»}u(p) (11) 
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Figure 1: One loop corrections to the elastic scattering of an electron by an external 
electromagnetic source 



where 



with 



rV,P) = ^iJd 4 kriSF(p'-k)YiSF(p-k)^iG aP (k) (12) 
P»tf,p) = j a iG a(3 (q)m^(q) (13) 

lV>p) = t^(p')tS F (p')Y + i^s F ( P )tj:(p) (14) 



ffl^(ff) = ^(-)Tif#kyPiS F (k + q)<fiS F (k) (16) 



<E(p) = ^ i jd i kiG a(3 (k) 1 a iS F (p-k) 1 ' 3 (15) 

(2tt)< 
and q = p' — p. 

We refer to eqs. (|12D, Cp~3|) , and flUf) as the Vertex, Polarization, and Leg con- 
tributions, which respectively correspond to diagrams (a), (b) and (c) plus (d). We 
also note that expressions ([12]) , fll5]), ( |16D represent the one loop corrections to the 
vertex, fermion and photon self energy parts respectively. 
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Given the form of the photon propagator it is convenient to introduce: 



1 



(17) 



where the subscript "0" denotes the (known) result coming from the standard part 
of the photon propagator, and for the part proportional to 7 2 in ([|) 



k 2 k 4 

In the remainder of this section we consider this part of the propagator only, omitting 
the label in the corresponding expressions. 

The procedure for evaluating the loop integrals is equivalent to that of standard (or 
metric) QED. We need to regularize them first and then renormalize the parameters, 
which include the THe\i parameter along with the fermion charge and mass. The 
regularization of the photon propagator is carried out using 

1_>_/ A2 dL 1 _ - 2 [* dL ( i 9 ) 

k 2 V {k 2 -L) 2 ' k 4 V {k 2 -Lf y } 

with the assumed limits // — > and A — > oo, and the parameter renormalization by 

the inclusion of the corresponding counterterms to each loop integral. Details about 

this procedure and the corresponding calculations are given in the appendix. We 

quote the final result for the loop integrals: 
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where we have implicitly assumed that ( pT|) is acting on free spinors. 
The Ward identity 

f^ = r>, P ). (23) 

is a consequence of gauge invariance, and therefore it holds even in the absence of 
Lorentz invariance. It is straightforward to check that ( p0| ) and (|2T| ) satisfy (|23|) . 

The evaluation of (|TT| ) is also straightforward once the loop integrals have been 
calculated. We just comment on the computation of the Leg correction, which is 
ambiguous since it contains terms like "0/0", which are indeterminate. To obtain an 
unambiguous result, we must explicitly introduce a damping factor, which is necessary 
for the correct definition of the initial and final "bare" states. Details of this adiabatic 
approach are presented in appendix B. The final result for the Leg correction is 



71 ( 



' ((3-p) 2 | 0-pP-g | l {(3-qf ~ 
i 2 m m 2 



m z 



+ P 2 



5 , , A . 9 133 1, ,m 

— ln(— Y H ln( — 

24 m 144 4 v 



(24) 



+ 
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6 m 



+ 



6 m 



Note that this part gives a contribution to the total amplitude that cannot be 
removed after renormalization. Furthermore, the gauge invariance of the Feynman 
amplitude which is manifest as 

q ■ A = (25) 

requires the presence of such terms, a condition that is not satisfied by the vertex 
contribution only. 

The final result for the scattering amplitude is 



(26) 
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(27) 
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p = (\LJL.E±JL £ + 1 )£\ (29) 

7r 1. 3 m. i fi m. 1 80 m. I m. \ 9A P> I m ) 



ir^3 m y6 m 

The various terms in Q2E|) distinguish the different contributions to the scattering 
amplitude. In (E7|) we group terms of order q 2 or higher. G M accounts for terms of 
order q at least, and 1^ for the gauge terms or those who give no contribution to the 
amplitude. Note that the remaining infrared divergence in can be understood in 
terms of soft photon radiation, analogous to the metric case. 

In the next section we will use the above results to compute the g — 2 anomaly. 

Ill (GM) g-2 

To lowest order the Feynman amplitude associated with the elastic scattering of an 
electron by a static external field is 

ieu{p') A(q)u(p) ■ (30) 

The radiative correction of order a to this process is given by 

ieu(p'){{l + £)Ao • A + 7 2 £A € • A}u(p) (31) 

where A represents the (known) metric result and represents the contribution 
from (p6|). 

In the nonrelativistic limit of slowly moving particles (|g| — * 0) and a static mag- 
netic field , it is straightforward to show that 

e fi('l) -» ~^B-S (32) 

eA r A -> eG ■ A (34) 

with C 1 given by (|2"8|), which is the dominant term as q — > 0. 

In order to simplify this contribution, we consider a constant magnetic field B, 
that is A = |r x B, in which case 

4 A^ -\^B -uo -u-B -ou 2 ) (35) 

where we have neglected the terms that mix the large and small spinor components. 
Similarly, we can show 

(3-q ^ A^ --(B-ua-u- B-ou 2 ) (36) 
2 
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and in the non relativistic limit 

iA ptZ~ iA p^_B.o (37) 
m 

and 

iA^-B-a . (38) 
If we put everything together in (f28|): 

eG ■ A—> — —-{B-ai— + — u 2 )--B -ua-u) (39) 
2m vr 42 12 ; 3 ' K J 

As a cross-check on the above result, we take the limit UiUj — > —5ij obtaining G-A —>■ 

—2Aq ■ A, which is the required limit consistent with the structure of eq. ( |18|) in that 

case. The previous result is the contribution of (|34]) to (|3ll) , which added to (|30|) , 

give us the relevant part of the Hamiltonian as 



H a = -{r S ■ B + Vj -uB-u} + 0{i 2 )0{a 2 ) = -r ij 'SiBj (40) 



with 



rs 2^ = ^ 2+ ? 1+ « 1 + ir (1+7s2))l} (41) 

where we have identified S = | , and u = u/|iZ|. The T parameters account for the 
coupling strength between the magnetic field and spin. We see that generalizes the 
gyromagnetic ratio of a fermion analogous to the manner in which the anomalous mass 
tensor generalizes the mass of a particle fl4l . We therefore identify the parameters 



r iJ = T5 l: > + r*-u l M J with the components of the anomalous gyromagnetic ratio tensor 
of the fermion in the class of THe/u, theories. 

Note that the presence of preferred frame effects induces a qualitatively new form 
of interaction between the spin and magnetic field which is quantified by r*. Here, 
instead of coupling with each other, they both couple independently to the fermion 
velocity relative to the preferred frame. This interaction stems purely from radiative 
corrections, and would be absent in any tree-level analysis of GMQED. 

Hence, eq. ( pEDp describes the interaction (as seen from the particle rest frame) 
between the particle spin and an external homogeneous magnetic field. From this we 
can extract the energy difference between electrons with opposite spin projection in 
the direction of the magnetic field as: 

eB r 



AE a -- 

2m 



g + g*u 2 cos 2 6 (43) 
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where B is the angle between the magnetic field and the preferred frame velocity. The 
influence of the radiative corrections (coming from g — 2 and g*) in (|43|) is negligible in 
comparison to the dominant factor of 2 in g. Since we want to single out the effects 
of the non-metric corrections, it is more interesting to study the precession of the 
spin or, more specifically, the oscillation of the longitudinal spin polarization. In the 
metric case, this frequency is proportional to the factor g — 2, and so it is a distinctive 
signature of radiative corrections. 

The observable quantity in the g — 2 experiments is actually the electron polariza- 
tion, which is proportional to the quantum mechanical expectation value of S, that 
is, (S). Using Ehrenfest's theorem, a quantum mechanical solution for the motion of 
(S) is obtained from the equation 

~^\r.f. = ~i[S, H a ] = Sx [TB' + r*(B' ■ u)u] (44) 

where the primed variables are referred explicitly to the particle rest frame (R.F.). 
Note that the preferred frame effect will show distinctly as a temporal variation of 
the spin component parallel to the magnetic field. 

In general we want to know the spin precession relative to some specific laboratory 
system, with respect to which the particle is moving with some velocity (3. This frame 
need not a-priori be the previously defined preferred frame, and so (3 ^ u. 

Since the THefi formalism does not change (locally) the fermion electromagnetic 
field interaction, we expect that a charged particle in the presence of an homogeneous 
magnetic field will satisfy the equation 

% = P x Q c (45) 
at 

with the cyclotron frequency f2 c = -^-B and 7 = (1 — (3 2 )~ l l 2 . Relating ( f4"4j ) to the 
laboratory system yields 

dS . dS 1 => -* , » 

~dt ^ ~ it ' + r x ^ m 

due to Thomas precession, with Qt = ~+i(df x This frequency is kinematic in 
origin and it is a consequence of the non-commutativity of the Lorentz transforma- 
tions. 

Relating the primed variables in fl4~4| ) to the laboratory ones by a Lorentz trans- 
formation gives 

dS, 
~dl 



Lab 



SxQ s (47) 
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with 

tt s = TB+ (1 -j)ti c + r*{B-u)u (48) 

where we have set E = and considered (for simplicity) the case of orbital motion 
perpendicular to the magnetic field (/3 • B — 0) in the above. Note that the spin 
precession about Q s is no longer parallel to the magnetic field (axial direction), but 
has a component parallel to u that comes from radiative and nonmetric effects. 

At this point it becomes necessary to define the preferred coordinate system. 
There are several candidates (such as the rest frame of the cosmic microwave back- 
ground) for this frame f2|. To study this issue it is sufficient to assume that the 
laboratory system (Earth) moves with a non-relativistic velocity (V) with respect to 
the preferred frame, and so we can identify 

u = V + (3 

In order to single out the effects of radiative corrections, we study the spin pre- 
cession relative to the rotational motion of the electron, that is: 

^-\ rot = SxQ D (49) 
at 

with tip = Sl s — i1 c and S = (S\_, S±, S\\), where the first two components are perpen- 
dicular to B (lower index) but parallel and perpendicular to (3 (upper index), and the 
last one parallel to B. In the following we refer to the difference frequency (Od) as 
the anomalous frequency (given its connection with the anomalous magnetic moment 
in the metric case). It is convenient to rewrite: 

n D = n a + n* a cos&(v ± + p) (50) 

with 

= ^ (g + 9*V 2 cos 2 9 - 2) B (51) 

and Q* = ^g*BV; where G represents the angle between V and the magnetic field, 
and V± the component of the velocity perpendicular to B. In Q a we group all the terms 
parallel to the magnetic field that contribute to the anomalous frequency (including 
nonmetric effects). The remaining terms perpendicular to B arise from nonmetric 
effects only, and produce a temporal variation of the spin component parallel to 
the magnetic field. This effect is absent in the metric case, and so represents a 
qualitatively new manifestion of possible EEP violation. 
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In general we are interested in solving ([49]) for the cases /3>>l / or/5<<l / so 
that 7(11) ~ 7(/3) or j(V), but is otherwise constant. Since Q* is proportional to £, 
we can perturbatively solve for each component in (|L9|) . Taking, for example, the 
initial condition 5(0) = S(3 we find 

S\_ = S cos Vt a t S± = S sin Q a t (52) 

Q,* D,* sin 29 

5|| = 5 , ^/3cos9(l-co S n a t) + 5' - + a Q [cos(^ a + Q c )t - 1] 

where we have chosen a coordinate system where B = z so that 

V - = B cos 9 + x sin 9, (3 = y cos Q c t — x sin Q c t (53) 

and assumed that any rotation related to 9 is negligible in comparison to other 
frequencies involved in the problem (fl a or Q c ). 

The fact that Q a was (in the metric case) proportional to g — 2, motived the very 
precise g — 2 experiments which were designed to specifically measure that anomalous 
frequency. We see that this frequency is modified by from its metric value by the 
additional terms present in fl5"T|). If we assume that the EEP- violating contributions 
to Q a are bounded by the current level of precision for anomalous magnetic moments 
||15|| , then the discrepancy between the best empirical and theoretical values for the 
electron yields the bounds 

|fe-| < 3.5 x KT 8 and |f e - - £ e+ | < 1(T 9 (54) 

the latter following from a comparison of positron and electron magnetic moments. 
For muons, a similar analysis yields 

|VI<10" 8 and lV-^+|<KT 8 . (55) 

Even though the accuracy of the muon anomaly is lower than the electron one, the 
slightly stronger bound in fl55|) arises because the experiments are carried out for 



high- velocity muons ||16|| . To our knowledge these bounds on violation of gravitational 
universality are the most stringent yet noted for leptonic matter. 

Newman et. al. analyzed the g — 2 experiments JTIj] in order to find new bounds 
for the validity of special relativity. They assumed that the parameter 7 involved in 
the electron motion had a different value (7) from that which arises kinematically (in 
Thomas precession and Lorentz transformations). The equivalent equation for (|5T|) 
is in that case 

n*™ = -(?-?! (56) 

m \2 7/ y ' 
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and by comparing with two electron g — 2 experiments, one at electron relativistic 
energy (/3 = 0.57) and the other nearly at rest ((3 = 5 x 1CT 5 ), they obtained the 
constraint £7/7 < 5.3 x 10 -9 . Our approach is qualitatively different from theirs, 
in that we assume 7 = 7 but include preferred frame effects in the evaluation of 
the anomalous magnetic moment. A similar analysis in our case yields the weaker 
bounds of |£ e | < 7 x 10~ 6 for electrons, and < 2 x 10~ 7 for muons. In the latter 
we used the g — 2 muon experiments carried at (3 — 0.9994 (7 = 29)[P!6[, and (3 = 0.92 
(7 = 12)0. 

Preferred effects not only modify the anomalous frequency according to fl5T|), but 
also induce oscillations in the spin component parallel to B. As stated above, this is 
a qualitatively new consequence of EEP violations due solely to radiative corrections 
in GMQED. Searching for such oscillations therefore provides a new null test of the 
EEP. We can estimate the magnitude of such effects by taking the temporal average 
of S\\ over the main oscillation given by Q a , which gives 

5= ^ ~ OTcose 7 2 (57) 

This effect is enhanced in highly relativistic situations, and can be estimated by 
considering a typical experiment with V ~ 10~ 3 . For electrons (3 ~ 0.5, and so 
5 e ~ 10~ n ; for muons (3 = 0.9994, yielding 5^ ~ 10 -8 . In both cases we used the 
corresponding present constraints for £ given above. 

The novelty of the S\\ oscillation suggests the possibility of putting tighter con- 
straints on the non-metric parameter, once appropriate experiments are carried out. 
The same goes for the analysis of Q a at different values of G (the angle between the 
magnetic field and the velocity of the laboratory system with respect to the preferred 
frame). The rotation of the Earth will turn this orientation dependence into a time- 
dependence of the anomalous magnetic moment, with a period related to that of the 
sidereal day. 

The previous analysis was concerned with effects related to spatial anisotropy. We 
turn now to considering possible violations of local position invariance. The position 
dependence in the former section was implicit in the redefinitions of charge, mass and 
fields. These quantities were rescaled in terms of the THefi functions, which were 
considered constant throughout the computation. LPI violating experiments are of 
two types. One of these entails the measurement of a given frequency at two different 
points in a gravitational field (where differences in the gravitational potential could be 
significant) within the same reference system. The other type involves a comparison 
of frequencies arising from two different forms of energy (i.e. two different clocks) 
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at the same point in a gravitational potential. We parameterize the gravitational 
dependence on a given frequency as: 



n = a 



1-U + EPUiA + ■■■ (51 



where C/y represents the external gravitational tensor, satisfying Ua = U, and the 
ellipsis represents higher order terms (going as either U 2 or velocity times U) in the 
gravitational potential or terms independent of it. 

The measured redshift parameter related to this frequency may be written as 

Z = AU(l-~), H = H«^ (59) 

where S* J will depend upon the specific frequency measured in the experiments. Note 
that this tensor is equivalent to the anomalous passive gravitational mass tensor 
introduced for the study of atomic transitions. 

In g — 2 experiments the relevant frequency is Q a , which describes the precession 
of the longitudinal polarization in the presence of a constant magnetic field. Using 
the THe/u, formalism (see eq.(|51~D) we obtain 



^ eB r , eB a 

tt a = — {g-2] + --- = - 

2m 2m ir 



-4 



+ • ■ ■ (60) 



where we have omitted terms proportional to velocities, which eventually will con- 
tribute as 0(v 2 U) terms at most. 

In order to carry out the loop calculation, the Tife/i dependence was absorbed 
into the definition of the parameters under the rescaling 

VI^-VL/cq m^-m^jTo/co a^a/e c (61) 

with Co = (To/Ho) 1 / 2 as the limiting speed of the massive particles, the subscript '0' 
denoting the Tife/i functions evaluated locally at X = 0. Although the product eB 
remains invariant under this rescaling, the expression for the constant magnetic field 
still depends on the THefi parameters once it is written solely in terms of atomic 
parameters. This can be seen clearly by considering the magnetic field produced by a 
long solenoid of length L, with N turns and carrying a current /. The gravitationally 
modified Maxwell equation to solve is: 

V x (fi^B) = An J (62) 

and so we find the non-vanishing magnetic field inside the solenoid to be B = 
AhhqIN / L. Again we assume that the THefi functions are constant throughout 
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the size of the experimental device. In terms of fundamental atomic parameters, L 
is proportional to an integer times the Bohr radius (the interatomic spacing), which 
is known to rescale as a — > aoeoCo/\/T^. If we now write I = J J ■ dS, where J 
can be expressed in terms of a density charge p in motion (v) through a volume V, 
and then relate the Bohr radius to each spatial dimension along with the limiting 
particle velocity Co to the velocity distribution v, we can show / — > Jv^o/ e o c o, and 
so B — > BjioTo/eQCQ. Along with ([61]) , this gives the position dependence of (^) to 
be 

° a = °^(^ (1 + 4 ) (63) 

with VL° a = eBa/2mii (recall £ = 1 — l//i eo c o)- 

Note that the THep functions are evaluated at some representative point of the 
system, which we have chosen to be the origin X = 0. In order to determine how Q a 
changes as the position of the system varies, we expand the THep functions in 



T(U) = T + Tfo ■ X + O(g ■ X) 2 (64) 

where go = Vt/|^ =0 , T = T\g =0 , and Tq = dT/dU\^ =0 . It is useful to redefine the 
gravitational potential U by 

U ^-\^9o-X (65) 
whose gradient yields the test-body acceleration g. This finally yields 

n a = a 



11 13 i 
1 - U + (— r - 4 A o)U] (66) 



6 u 6 



where we have rescaled again according to ( |61]) , and omitted terms proportional to 
£, since the main position dependence parameterization is given in terms of: 

r _ 2T 0/ e d , T H N A _ 2T o,Mo , T H x 

By comparing eq. (|66|) with ( |58f) , we can identify 

5 s - 2 = ^r - ^A (68) 
6 6 

as the LPI-violating parameter. Note that this depends on the anomalous frequency 
related to the longitudinal polarization of the beam. It is also species-dependent, 
with the value of To and Ao for the electron differing from that of the muon. A search 
for possible position dependence of anomalous spin precession frequencies provides 
another qualitatively new test of LPI sensitive to radiative corrections. 
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Actually the most precise g — 2 experiments for electron measure the ratio a = 
Q a /Q c at non relativistic electron energies {(3 ~ 10~ 5 ), and so Q c ~ eB/m. This is 
interesting because by following the former parameterization we can write: 



tt c = {l° c \l - U + (2T - A )U 



or by taking the ratio of (|66|) to (69): 

a = a°(l + UE a ) 



7i T o + l A o 
6 6 



(69) 



(70) 



and then by identifying a with the most precise experimental value |H| and a with the 
theoretical one [|T5|1 , we can constrain through the resulting theoretical/experimental 
errors \UE a \ < 3 x 10~ 8 . This result is sensitive to the absolute value of the total 
local gravitational potential M, whose magnitude has recently been estimated to be 
as large as 3 x 10~ 5 due to the local supercluster |Q. Hence measurements of this 
type can provide us with empirical information sensitive to radiative corrections that 
constrains the allowed regions of (ro,Ao) parameter space, giving in this case: 



l^r + Ja | 

6 6 



< 10 



-3 



(71) 

and so nothing conclusive is 



For muons the analogous constraint is \UE®\ < 10~ 5 , 
obtained. 

We pause to compare this result to an analogous result obtained for hyperfme tran- 
sitions (maser clocks). In this case the baryonic and leptonic gravitational parameters 
appear simultaneously. The atomic hyperfme splitting comes from the interaction be- 
tween the magnetic moments of the electron and proton (nucleus). The proton metric 
appears only in the latter, and so it does not affect the principal and fine structure 
atomic energy levels. Non- metric effects imply a shift in the hyperfme energy Ehf 
which is @ 

AE hf = E hf {\ - U B ) + £ hf U B E hf (72) 

with 

E hf = 3V B -A B + A (73) 

where Ub,T b and A B are the baryonic analogues of the parameters appearing in d67|). 
In (|72|) we rescaled the atomic parameters to absorb the THefi functions and chose 
units such that c B = 1. The quantity A is given by 



A 



2^ 

rpl 



2 (Ek _ ^o- 



H 



B 



T B T 



(74) 
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and would vanish under the assumption that the leptonic and baryonic THe^i parameters 
were the same. The gravity probe A experiment employing hydrogen-maser 



clocks, was able to constrain the corresponding LPI violating parameter related to 
hyperfine transitions, obtaining 

\E Hf \ = \3T B -A B + A\ <2 x 1(T 4 (75) 

for the most stringent bound to date on "E H * . 

We note that a similar analysis could be carried out for the energy shift defined in 
(f43|), which can be used as a frequency test to look for position or frame dependence. 
This can be done by following the same procedure as for atomic energy shifts, where 
the anomalous passive and inertial gravitational tensor are introduced in order to 
relate non-metric effects to redshift and time dilation parameters. Since radiative 
corrections are irrelevant in that energy shift, we omit that procedure here. 



IV Concluding Remarks 

Refined measurements of anomalous magnetic moments can provide an interesting 
new arena for investigating the validity of the EEP in physical systems where radiative 
corrections are important. We have considered this possibility explicitly for the class 
of non-metric theories described by the THen formalism. The non-universal character 
of the gravitational couplings in such theories affects the one loop corrections to the 
scattering amplitude of a free fermion in an external electromagnetic field in a rather 
complicated way, giving rise to several novel effects. 

An evaluation of the one-loop diagrams reveals that the leg corrections, which 
in the metric case give no contribution to the total amplitude after a proper renor- 
malization of mass and spinor field, provide contributions which cannot be removed 
after renormalization. Moreover they are essential in ensuring the gauge invariance 
of the scattering amplitude, which is not fulfilled by the vertex correction alone. The 
consistency of the calculation is verified explicitly through the Ward identity, which 
furnishes a cross-check between the fermion self energy and the vertex correction. The 
non metric corrections to the scattering amplitude also have an infrared divergence, 
which could be understood in terms of inelastic soft photon radiation, as in the metric 
case. This does not affect the term associated with the anomalous magnetic moment. 

The presence of preferred frame effects induces a new type of coupling between 
the magnetic field and the spin as described by (^0[). This interaction stems purely 
from radiative corrections, and generalizes the gyromagnetic ratio of a fermion to a 
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tensorial coupling described by T^. We emphasize that qualitatively new information 
on the validity of the EEP will be obtained by setting new empirical bounds on this 
coupling, as it is associated with purely leptonic matter. Comparatively little is known 
about such empirical limits on EEP- violation relative to the baryonic sector j2J, for 
which previous experiments have set the limit ||] \£g\ = \1 — c 2 B \ < 6 x 10~ 21 where 
cb is the ratio of the limiting speed of baryonic matter to the speed of light. We can 
therefore safely neglect any putative effects of £b in our analysis. 

Consequently, discussion of a g — 2 contribution to the magnetic moment no 
longer makes sense, and we instead refer to the anomalous frequency as the main 
connection with experiment. Note that this frequency, defined as the relative electron 
spin precession with respect to its velocity, comes from radiative corrections and it 
becomes proportional to g — 2 in the metric case. This frequency shows an explicit 
dependence on both the preferred frame velocity and its relative direction with respect 
to the external magnetic field. There is also a dependence on the electron velocity, 
which makes the other contributions negligible at relativistic electron energies. Two 
g — 2 experiments on the electron (one at relativistic energies and the other almost 
at rest) may then be used to limit the preferred frame parameter to be no larger 
than 10 -5 , analogous to the work of Newman et ai. Constraining any possible EEP 
violation to be no larger than the present discrepancy between theory and experiment 
we found the most stringent bounds for £ yet obtained for leptonic matter, as given 
in (H) and fl55|). 

We expect that new experiments which probe the anisotropic character (or angular 
dependence) of the frequency could be used to impose stronger limits in different 
physical regimes. For example, as the Earth rotates, the spatial orientation of the 
magnetic field changes - this should in turn diminish the experimental errors involved 
in the comparison between two energetically different g — 2 experiments. 

The relativistic generalization of the spin polarization equation (f47|) , followed the 
same procedure as for the metric case, where non-metric effects where included in 
the interaction only (eq. (flOP). This yields an equation of motion for the spin (as 
seen from the rest frame) which is qualitatively different from that expected from 
its classical counterpart, where the angular momentum rate is related to the torque 
applied on the system. This approach for dealing with violations of Lorentz invari- 
ance is dynamical; from a kinematical viewpoint we assume that standard Lorentz 
transformations relate coordinates and fields from one system to another. 

Perhaps the most remarkable feature of the non-metric effects is that of the os- 
cillations of the component of spin polarization parallel to the magnetic field. Since 
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this component remains constant in the metric case, an experiment which searches 
for such oscillations is a new null test of the equivalence principle that is uniquely 
sensitive to radiative corrections in the leptonic sector. Hence an empirical investiga- 
tion of its behavior will provide qualitatively new information about the validity of 
EEP, and could constrain even further the limits on the preferred frame parameters. 

Finally, we analyzed the behavior of the anomalous frequency in the context of 
redshift experiments, which can put constraints on the LPI-violating parameters (To, 
A ) once the corresponding experiments are carried out. This region of parameter 
space is qualitatively different from that probed by either Lamb-shift or hyperfme 
effects. In the electron sector a bound on the magnitude of lfE a can be obtained by 
demanding that it be no larger than the error bounds in the discrepancy between 
the experimental and theoretical values of the ratio a = Q a /Q c . Assuming the local 
potential to be as large as that estimated from the local supercluster, we obtain a 
bound on |S a | that is comparable the limit on an analogous quantity in the baryonic 
sector obtained from redshift experiments However this latter experiment is 



proportional to changes in the local potential, which are ~ 10 -10 . More direct limits 
on | S a | must be set by performing a similar sort of redshift experiment on anomalous 
magnetic moments [pH[] . The logistics and higher precision demanded by such an 
experiment will be a major challenge to undertake. 
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Appendices 
A Loop integrations 

We show the main steps leading to Eqs. flU), ©, and (H). Details are given 
throughout the computation by considering only the first term of the photon propa- 
gator (|18|), that is 

Gl = ^ + --- (Al) 



with the remaining term in (|18|) contributing in a similar manner. 
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We solve for the fermion self energy by replacing (XI) in (I5[), and using ([I 
along with the Feynman parameters 



dz z 



a?b Jo [az + b(l — z)f 



After integrating we obtain 



S(p) = / 0(jfe + m) /Hn 

47T J 



'm 2 (l - z) 2 + zA 2 - Az(l - z) 



(A2) 



m 2 (l — z) 2 — Az(l — z) + z/i 2 y 

with A = p 2 — m 2 . We consider A/m 2 << 1, and expand the above to obtain after 
some manipulation 



(A3) 



— m 



3Vl_ ln( A. 



2 V2 



+ 2^|^(^ - m) f ln(-) 2 - 3] } + 0((1- m) 2 ) + 



where we have kept the leading terms as \i — > and A — ► 00, and O 
for the terms satisfying 

S F (p)0((p'-m) 2 )u{p) = 0. 
We renormalize S(p) by subtracting 

5E = 5m + 8£ 0/3 ■ p 



m) 2 ) stands 



(A4) 



where the counterterms respectively account for mass and THefi -parameter renor- 
malization. 

Choosing the counterterms so that 



and so 



5m 

the regularized result is then 
E(p) = 

7T 



u(p)T,(p)u(p) = 
a (3 2 (I 



va- 



ra 



4vr 2 V2 
4^ V2 + W 



I ln( - )2 _3 
2 > 2 



m 



(A5) 
(A6) 



/? 2 



1, ,A. 2 1 
8 ra lb 



+ 0((^-m) 2 ) 



(A7) 
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Note that the remaining ultraviolet divergence related to this term could be removed 
after charge renormalization. We find it convenient to leave it in order to cross-check 
the calculation, since a similar term from the vertex part should cancel it, thereby 
removing the divergence from the resulting scattering amplitude. 

The evaluation of the vertex function follows a similar procedure, giving the result 



a r dx ( 
~ / ~V 

IT J P X V 



^•^•p'(ln^-2) + ^-ln^ 



x\ (P ■ VI 11 ft + P ■ p' Pi 1 ") ~ (P ■ Pf V + P ■ P' h") (1 - x£ (A8) 



+ fife* ■ PPx ~ \^Pl(\ ~ ln ^) " (! - X )P ■ PP" ~ *P' -PP") 



+ 



■ T ' 1X 1 rxK 2 
with p x = xp' + (1 — x)p. Since ( |X8| ) is acting on a free spinor, we can use 

2 2 (^ \ 2 

p x = m — x(l — x)q , 

with q = p' — p, and so expand 



— = l + x(l-x)^ + 0(q i ) 



p: 



which after some algebra reduces ( |A8| ) to 

7T <- 



16 8 m 



+ 



,m 



P 2 



m 2 



P-pP-qfl 



,m N 



+ P 
m 



16 



4 m 8 m 



m m 

2 



+ 



5 q 
24 m 1 



2 m m 8 m 

p 2 q»- 



lP-p Iff-gW 
Am 3 m I m 



(A9) 



-} + 0(q 3 ) 

where the vertex function has been renormalized by subtracting a term like 

5T^ = 5£ W (A10) 

with <5£o is given by ( |A6|) . We recall that gauge invariance forces this coefficient to 
be equal to the one participating in the renormalization of the fermion self energy. 
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B Adiabatic hypothesis 



In order to describe how self energy effects convert the incident electron from a bare 
particle to a physical one, it is convenient to introduce a damping function, g(t), 
which adiabatically switches off the coupling between fields, such that the interaction 
lagrangian is replaced by 

d = eg{t)^{x) A(x)i/;(x) (BI) 

It is assumed that the time T over which g(t) varies is very long compared to the 
duration of the scattering process. In momentum space 

g (t) = J G(n ) exp(ifi • x)dtt (B2) 

with Q = (fio> 0), and g(0) = 1. It is supposed that G(Qq) is almost a delta function, 
being large for Q in a range of about T _1 

In the presence of an external field A^, eq. fll4|) will now read 

l-a—> J G(n )G(%)dn d% A(p-p-n-n')tS F ( P -n-n')ii:(p-n) + --- (B3) 

where ■ • ■ represents the equivalent second term from 

As T — > oo, and Qo, - * 0, the fermion propagator reduces to 

S F = l - ~ ^ + m (B4) 

where we used p 2 = m 2 . This implies that we can expand S up to order Q only, since 
higher terms vanish after taking the previous limit. Here we employ the relation 

1 \ + \bU- (B5) 



A-B AAA 
to expand 

E(p - O) ~ E(p) - ^-Jy • (B6) 
After renormalization, £ takes the form 

Y,{p) = {p'-m){A + B{(3-p) 2 ) + C{(3-p){ m) + (jf- m) ^} + • ■ • (B7) 



where the constants A, B, and C can be obtained from eq. (|20|j . 

Let us introduce = O + Q', and symmetrize fl by ~(Q + f2') in (|B6|) , to write 

s f (p - a - n')E(p - n) = ^t^(s(p) - m 
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which after using ( [B7| ) can be written as 

l(A + B((3-p) 2 ) + Cp-p0 (B9) 

where we have used that S(p — Q) is acting on a free spinor, and therefore terms of 
the form (jrf— m)u(p) vanish. Now, the final evaluation of (|T4]) follows directly from 
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